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The coupling of two s-wave superconductors through a small magnetic dot is discussed. Assuming that the 
dot charging energy is small compared to the superconducting gap, E c <C A, and that the moment of the dot is 
classical, we develop a simple theory of transport through the dot. The presence of the magnetic dot will position 
Andreev bound states within the superconducting gap at energies tunable with the magnetic properties of the 
dot. Studying the Josephson coupling it is shown that the constructed junction can be tuned from a "0" to a 
"7r" -junction via a degenerate two-level state either by changing the magnetic moment of the dot or by changing 
temperature. Furthermore, it is shown that details of the magnetic dot can be extracted from the sub-harmonic 
structure in the current-voltage characteristics of the junction. 
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If a superconductor is exposed to magnetically 
active impurities 0,| or materials § the super- 
conducting state is modified. Coupling two super- 
conductors through a magnetically active barrier 
may lead to what is known as a "7r" -junction 0, 
a junction for which the ground state has an in- 
ternal phase shift of ir between the superconduc- 
tors across the barrier. If the barrier is extended 
to an S/F/S-structure, a ferromagnetic (F) layer 
sandwiched between two superconductors (S), the 
critical current will oscillate as the thickness of F 
is varied || . Recently, this behavior was reported 
in a mesoscopic S/F/S-structure and a system- 
atic change from a "0" to "it" Josephson junction 
as a function of the F-layer thickness was seen 
g . In this paper we present a different scenario 
which gives rise to qualitatively the same physics, 
namely, two s-wave superconductors connected 
through a small magnetic dot. In this case the 
effect of the magnetic dot arises from the spin- 
active contact at the interface between the dot 
and the superconductor in contrast to the inter- 
nal exchange field of the fcrromagnet as in Refs. 

Hi- 

In Figure |l| we show a cartoon of our sys- 
tem: two conventional superconductors brought 




Figure 1. A schematic illustration of the super- 
conductor/magnetic dot/superconductor contact 
considered. The size of the dot, d is small com- 
pared to the superconducting coherence length £o- 
It is also assumed that the moment, /J, is large, 
i.e. in the classical limit. 

in electrical contact through a magnetic grain. 
The size of the grain, d, is small compared to 
the superconducting coherence length, £o, so that 
<i/£o "C 1 holds. The dot is characterized by a 
large classical magnetic moment, /2, and despite 
the small size of the grain we shall assume that 
charging effects on the dot can be neglected, i.e. 
Ec < A. In our approach the effect of the mag- 
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netic dot enters via boundary conditions for the 
quasiclassical Green's function. In this spirit the 
contact over the magnetic dot is described by the 
following S-matrix 

S=(; _ r *)exp(i|^) (1) 

were r and t are the usual reflection and trans- 
mission coefficients. The spin-dependence enters 
via a rotation of the quasiparticlc spin by an an- 
gle O around the moment axis, a^. The angle 
O is a phenomenological parameter that quanti- 
fies the degree of "spin-mixing" induced by the 
moment p, and varies between (no spin-mixing) 
and 7r (strong spin-mixing) . The calculations 
presented here are done within the quasiclassical 
theory as presented in Refs. [||J^] and references 
therein. 

Local quasiparticle density of states. 

One of the most pronounced effects of the mo- 
ment, /2, on the superconducting density of states 
(DOS) is the formation of Andreev states at 
£b(0) = ±Acos(0/2). The two signs refer to the 
two " spin bands" on which the pairing amplitudes 
are ~ (ipfipi) (+) and ~ (VmV't) ( — ) respectively. 
The formation of these Andreev states is quite 
easily understood since the moment will give re- 
flected quasiparticles a phase shift of exp(±i0). 
The sign of the phase shift depends on the spin of 
the quasiparticle: (+) for spin up and (— ) for spin 
down. Introducing the moment axis, oyx, breaks 
the rotational symmetry in spin-space of the su- 
perconductors and pins down the direction of the 
spin-quantization axis to be parallel to &fj,. 

Formation of Andreev states is a signature of a 
reduced superconducting order by pair-breaking. 
The Andreev states are spatially localized at the 
magnetic moment and their spectral weight de- 
cays exponentially in to the superconductors over 
a distance given by £(eb) = vp/2y/ A 2 — e 2 b , 
where vf is the Fermi velocity in the supercon- 
ductor. This gives that superconductivity heals 
to its bulk value over a distance ~ S,(sb) into the 
superconductor. However, in the limit d/^o < 1 
the pair-breaking effect of the dot on the super- 
conducting order parameter is small and we can, 
to good approximation, assume a constant or- 



der parameter up to the contact. Another con- 
sequence of the relative smallness (d/£o ^ 1) of 
the dot is that the broadening of the Andreev 
state by elastic impurity scattering is suppressed. 




Figure 2. The density of states (DOS) at a mag- 
netic dot for different values of in (a). The 
bound states appear at eb = ±Acos(0/2), the 
different signs refer to two the spin bands (see 
text). Here we plot eb = +Acos(0/2). In panel 
(b) we display the DOS for the momenta ±p, on 
the same spin-band as in (a), with = 37r/4 for 
various values of phase difference x over the junc- 
tion, V = 0.9. 

It is important to note that the positions of 
the bound states do not disperse with momen- 
tum p but do only depend on the angle 0. This 
makes the position of states Eb robust and the 
DOS shown in Figure ||.a are angle averaged ones 
at the dot. Also important to note is that the 
DOS in both superconductors are symmetric and 
Eb coincide in both sides. This makes the DOS 
independent of the transmission probability, T>, 
as long as there is no phase difference \ applied 
across the junction. 

Tunable Josephson-junction properties. 

It is well known that the Josephson current is car- 
ried by junction Andreev states p|-|Ic|]. In case of 
a magnetic dot the junction Andreev states split 
into four states, labelled by their momentum di- 
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rections ±p and by spin-band, as a phase differ- 
ence x is applied across the junction 0. We have 



e B ( X ) = ±A 



cos 2 (f)-2?cos(e)sin 2 (f) 



±VPsin(9)sin(|) J I -X>sin 2 (f) 



(2) 



which now not only depends on but also on \ 
and T>. In Figure[2|b we show the DOS for the two 
opposite momentum directions, +p and — p, on 
the positive spin-band with = 37r/4, for various 
X at T> = 0.9. As seen in Figure |].b the Andreev 
states in the DOS splits equally for both +p and 
p with applied phase difference x as given by 
Eq.(H). The spectral weight of the bound states, 
w(sb(x, P))i however, is not symmetric in mo- 
menta, i.e. w(s B {x, -p)) ^ w(e B (x,+P))- This 
leads to, after summing up all contributions, a 
net Josephson current that may have additional 
zeros in the phase interval and tt. 

The sensitivity of the quasiparticle spectra in 
Eq.(||) to the degree of spin-mixing, O, and to 
the transmission probability T> at finite phase dif- 
ference, makes the magnetic dot a very tunable 
Josephson contact between the two superconduc- 
tors it connects. In d-wave superconductors, as 
the high-T c cuprates, Andreev states at zero- 
energy lead to anomalous Josephson properties 
pTp^ ]. In an ideal junction between two d-wave 
superconductors the zero-energy states will cause 
the critical current to diverge ~ 1 /T at low tem- 
perature (TJ]. In addition, the energetics of the 
junction may change from being a "0" -junction 
at temperatures just below T c to a "7r" -junction 
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at low temperatures T <^T, 

As demonstrated in Refs. |6|,7| the Josephson 
coupling of two s-wave superconductors through 
a magnetically active barrier shows much of the 
same anomalies as seen in the d-wave / d-wave con- 
tact. In Figure (3|a we show the critical current 
as a function of temperature for several values of 
O ranging for to tt. At low temperatures and 
for O close to tt the critical current is roughly an 
order of magnitude larger than the value of a non- 
magnetic dot, i.e. at O = 0. In fact, it is a finite 
transparency, in Figure ||.a put to I? = 0.1, or 
an inelastic or phase-breaking life time, Ti ne i )P b, 
depending on which gives the largest broadening 
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Figure 3. The critical current vs temperature of a 
magnetic-dot Josephson contact is shown in (a). 
The spin- mixing angle, O is varied from to 7r 
(top to bottom) in steps of 7r/20. The trans- 
mission probability V is 0.1 (T> — |i| 2 ). In pan- 
els (b) and (c) we show the current-phase rela- 
tion (b) and the junction energy vs. phase (c) 
for O = 37r/4 at temperatures just around the 
switching temperature T sw ~ 0.12T C , highlighted 
by the dashed circle. At T sw the junction changes 
state from a "tt" to a "0" -junction. As seen, 
the junction goes through a degenerate 7r-periodic 
state at the switching temperature. 

of the junction states, that cuts of the growing 
critical current as T — > (|J^] . 

In Figures ^.b-c we show the current-phase and 
the energy-phase relations of a junction with O = 
37r/4 and T> = 0.1. The temperature is swept over 
a small interval including T sw . At T = T sw the 
junction changes its energy state from being a 
"0" -junction at low temperatures to being a "ix"- 
junction at T > T sw . This transition goes over 
7r-periodic current-phase relation at T — T sw and 
here the "0" and "tt" states of the junction are 
degenerate. A similar transition may be achieved 
at fixed T and instead sweeping O. 

Current-voltage characteristics. 

As introduced, the spin-mixing angle is a phe- 
nomcnological parameter in the same spirit as the 
junction transmission probability is. In principle 
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Figure 4. The current-voltage characteristics for 
the magnetic Josephson contact shown for © = 
0, 7r/2, 37r/4, and 6 = n in panels (a) through 
(d). In each panel we show the IV-curve for 
different transparencies as labeled in panel (a). 
The IV-curves in panel (a) correspond to an or- 
dinary single-channel point contact between two 
s-wave superconductors. The appearance of An- 
dreev states at eb — ±Acos(0/2) changes the 
sub-gap structure of the magnetic dot (b-d) as 
described in Eq.(|^) and here marked by the ver- 
tical dashed lines. In each case the current is 
given in units (GatA/c) where Gm is the normal 
conductance. 



one should be able to give a microscopic justifi- 
cation of the value of 9 by calculation along the 
lines in Ref. ||. In a practical realization of a 
magnetic dot it is more likely that an effective O 
is measured and it is important when characteriz- 
ing the junction to extract the value of O by com- 
paring measurements with calculations. Pinning 
down and T> can be done by measuring several 
junction specific properties such as the Joseph- 
son current, the current-voltage characteristics 
and the noise spectrum of the junction as was 
recently done in Al point contacts jL3|. Here we 
compute one of these properties, the dc-current- 
voltage curves, for a single channel contact as a 
function of transparency, T>, and spin-mixing, 0, 
using the methods described in Ref. 0. 

At a constant voltage bias, \eV\ < 2A, applied 



over a junction between two superconductors the 
only way a dc-current can flow through the con- 
tact is via multiple Andreev reflection (MAR) 
|l4|| . This simply means that in order to transfer 
a quasiparticle at energy e = — A in the super- 
conductor at bias — \eV\/2 to energy e = A in 
the superconductor at bias +\eV\/2 the quasipar- 
ticle must undergo 2A/|eV r | Andreev reflections. 
At each reflection event the quasiparticle is ac- 
celerated and shifted up in energy by an amount 
\eV\. In the current- voltage characteristics this 
shows up in an abrupt increase in the current 
each time the voltage is swept over a value where 
mod(2A, \eV\) changes. This holds true for the 
conventional superconductor for which there are 
no states within the superconducting gap [ p~5| , p~6| . 
The resulting sub-gap structure in the current- 
voltage curve for = is shown in Figure^. a for 
several values of T>. 

Introducing Andreev states within the gap, by 
a finite value of 0, we modify the MAR-picture 
above. Now, having the bound state es(0) = 
±Acos(0/2) and starting at a large voltage bias 
|eV|~2A, the first voltage below 2A to connect 
two points in the quasiparticle spectra, on ei- 
ther side of the junction, with a finite DOS, and 
thus to give a contribution to the current, is at 
\eV\x = (1 + cos(0/2))A. There are two pos- 
sible processes giving current at \eV\\ per spin- 
band. For the band with e B (Q) = +Acos(0/2) 
we have: (i) connecting an energy at —A with 
£b(0), i.e. an electron- like quasiparticle tun- 
neling into the bound state from the continuum 
states below — A. (ii) connecting £_b(0) with the 
continuum states at —A, i.e a hole-like quasipar- 
ticle tunneling from the bound state into the con- 
tinuum states below —A. For the spin-band with 
£s(0) = — Acos(0/2) the same processes occur 
with the modification that the quasiparticle (hole) 
start in the bound state and ends up in the contin- 
uum at +A. The next "spectral current channel" 
to open is at voltages \eV\2 = (1 + cos(0/2)) — 
which, via a single Andreev reflection off the op- 
posing superconductor, connects finite DOS in 
the same electrode. Adding more and more spec- 
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tral current channels, |eV| n , we arrive at 
|eV| B =(l + «»(§))£ 



(3) 



which gives the location of the strong features 
in the current-voltage curves at different spin- 
mixing angles, O. This is clearly demonstrated in 
Figure fl| where the position of the sub-gap struc- 
tures described by Eq.(||) are marked. 

Why do you not see structure in \eV\ n = 2A/n 
as well? In principle, one could think that the 
usual subharmonic gap structure should be su- 
perimposed on this new one. The answer is that 
you still have the usual processes between — A and 
+A, but now the DOS at the gap edges is rather 
small. Thus, opening the possibility for one of 
these processes, its probability is small (although 
its contribution will increase with voltage). 

Finally we should mention that the IV-curves 
in Figure [| are calculated with a very small intrin- 
sic broadening, l/2ri„ e ; :P b = 1CP 4 A. Increasing 
l/2Ti ne i t pb smears out the sharp features at low 
voltages in to a large current peak. This is par- 
ticularly strong when sb ~ 0, i.e. for O close to 
7r. At = 7r the IV-curves are similar to those of 
a symmetric, 45-degree rotated, d-wave/d-wave 
junction 17 7j] ■ 



Conclusions and acknowledgements. 

In conclusion we have presented a simple the- 
ory to handle transport through magnetic dots 
connecting two conventional superconductors in 
the limit where charging effects may be neglected. 
We show that the Andreev bound state spectra 
is tunable with the properties of the dot. This 
allows realizations of Josephson junctions that 
are either "0" or "^"-junctions. Furthermore, 
we look at ways to extract information of the 
dot from experimental data and in detail study 
the current-voltage characteristics. The presence 
of Andreev states in the DOS allows for current 
through resonant tunneling flg}| and modifies the 
sub-gap structure in the IV-curves compared to 
the case of a conventional SIS-contact. 
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was founded by the Swedish Research Council 
(VR) (M. A. and M. F.) and by the DFG project 
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